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^ ■ Abstract 

00 , 

•/^ , We extend our analysis of transverse single spin asymmetry (SSA) in electroproduction of J/V' 

■»H- ' to include the effect of the scale evolution of the transverse momentum dependent (TMD) parton 

o' 

en ' distribution functions and gluon Sivers function. We estimate SSA for JLab, HERMES, COMPASS 

and eRHIC energies using color evaporation model of charmonium production, using an analytically 

S^ I obtained approximate solution of TMD evolution equations discussed in literature. We find that 

.^.' there is a reduction in the asymmetry compared with our predictions for the earlier case considered 

by us, wherein the Q^ - dependence came only from DGLAP evolution of the unpolarised gluon 

densities and which used a different paameterisation of the TMD Sivers function. 
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I. INTRODUCTION 

The factorization theorems of Quantum Chromodynamics(QCD), which enable us to 
calculate the cross sections of hadronic processes as a convolution of partonic cross section 
with parton distribution functions(PDFs) and fragmentation functions(FFs), are based on 
a collinear approximation. In this approximation, the intrinsic transverse momentum is 
integrated over in the definition of PDFs and FFs. However, this collinear factorization at 
leading twist is unable to account for the Single Spin Asymmetries (SSA's) that have been 
observed in processes involving scattering off polarized hadrons \Mi^- One of the approaches 
that have been used to explain these SSA's is based on the TMD factorization formalism 
wherein the PDFs and FFs depend on intrinsic transverse momentum k± also in addition 



to the momentum fraction variable x 



Ml- 



TMD factorization is particularly useful for 
describing the processes that are sensitive to parton's intrinsic transverse momentum. A 
formalism for TMD factorization has been developed by Collins, Soper and Sterman [9|] 
and has been used to study Drell-Yan (DY), Semi inclusive deep inelastic scattering(SIDIS) 
and back-to-back hadron production in e"'"e~ annihilation at small transverse momentum 
[lO(]. There has been huge interest amongst both theoreticians and experimentalists in 
the study of the TMD PDFs, TMD FFs and the resulting SSA's with the aim of probing 
the transverse momentum distribution of partons inside hadrons and the spin structure of 
hadrons. Henceforth, we will refer to TMD PDFs and TMD FFs collectively as TMDs. 

One of the important TMDs is the Sivers function which describes the transverse momen- 
tum distribution of an unpolarized parton inside a transversely polarized hadron. Estimates 
of Sivers asymmetry have been given in Drell Yan (DY) process. Semi-inclusive Deep In- 
elastic Scattering (SIDIS) and D-meson production based on a generalized factorization 
formula involving TMD PDFs and TMD FFs. In a previous work [ll], we proposed Jf-ip 
production as a probe of gluon Sivers function and gave estimates for SSA's in low virtuality 
electroproduction of J/ip at JLab, HERMES, COMPASS and eRHIC energies using color 
evaporation model of charmonium production. Investigation of SSA's in charmonium pro- 
duction as a possible tool towards resolution of charmonium production mechanism puzzle 
has been discussed within the NRQCD framework some time back [12]. It was shown in 



Ref. I4], using general arguments, that SSA in heavy quarkonium production is sensitive 



to the color configuration of the qq pair produced at short distance. Based on an analy- 



sis of initial and final state interactions , the authors argued that in NRQCD framework, 
asymmetry is non-zero in e-p collisions only in color octet model and in p-p collisions it 
is non-zero only in color singlet model. In our earlier work in Ref. ll| we performed a 



phenomenological study of the SSA in charmonium production. We predicted non- zero 
asymmetry in the process e -f p^ — ;■ e + J /ip + X using color evaporation model (CEM) of 
charmonium production [l3|. In the older work, we had used parameters of Sivers function 
fitted to SIDIS data in Ref. [ij]. Earlier parameterizations of TMDs like the ones we had 



used 



111 ] had been obtained by assuming a parton model picture of TMD factorization and 



directly fitting IJ, |l5| the calculated cross sections to experimental data from HERMES 
and COMPASS. However, the Q^ - evolution of TMDs had not been taken into account 
in obtaining these parameterizations. Our earlier work which used these, thus could not 
include effect of evolution of TMDs. 



In present work, we improve our previous estimates by taking into account the Q^ evo- 
lution of TMD PDFs and the Sivers function and energy dependence of Sivers asymmetry 
has attracted a lot of attention in recent past 16l-l2l|. A TMD factorization framework 
taking into account Q^ evolution of TMDs has been proposed [l6|. This framework has 
been used to obtain evolved TMDs from fixed scale fits. A simple strategy for extracting 
TMD evolved Sivers function has been proposed and a best fit of the SIDIS asymmetries 
has been performed taking into account Q^ -dependence of TMDs 22j using data from 



Ref. [23|]. The authors have also compared these fits with those obtained earlier [l4j. In 
fact, the TMDs evolve much faster than the integrated PDFs, for which the scale evolution 
is given by the DGLAP equation. In this work, we apply the formalism of Ref. 22| to 
J /ip electroproduction to improve our earlier estimates, taking into account the effect of Q^ 
evolved TMDs. In the current investigation, we restrict ourselves to the use of analytical 
formulae for the evolved TMDs, given in Ref. [22|. 



In Section II, we review the formalism 



or J JTp production in CEM and also summarize 
the TMD evolution formalism of Ref. 22|. In Section HI, we present our estimates for 
asymmetry obtained by taking into account the TMD evolution and scale dependence of 
Gaussian width of PDFs. Section IV contains summary and discussion of our results. 



II. FORMALISM 



A. ASYMMETRY IN J/^ PRODUCTION 



We consider the leading order (LO) parton model cross-section for low virtuality electro- 
production (photoproduction) of J /ip within color evaporation model. According to CEM, 
the cross section for electroproduction of J /ip using WW approximation is given by, 



o 









where rric is the charm quark mass and 2mu is the DD threshold. 

To calculate SSA in scattering of electrons off a polarized proton target, we assume a gen- 
eralization of CEM expression by taking into account the transverse momentum dependence 



of the Weizsacker- Williams function and gluon distribution function [11 1 



Q^fCC 



= / dXy dXg [d k±yd k±g\ jg/pt{Xg,k±gjj^/f.[x^,k±yj—-—— [2) 

where M^ = M^g. The difference in da"^ and da^ is parameterized in terms of the gluon 
Sivers function 

da^ — da^ = / dx^ dxg d'^k±^ d^k±g A'^ fg/pt{xg, k±g) f^/dxy, k±y) da'^^^'^'^ (3) 

where da is the elementary cross section for the process 'jg — )■ cc given by 

1 d^Pcd^Pc 1 r4/ \ 7T7 [2 /^\ 

^^ = 77^ 77^77^ 777^N2 ^ (^7 + Pg ' Pc ' Pc) \ M^g^^c \ ■ (4) 



2s 2Ec 2Eg,_ f27r' 



Following the procedure in Ref. 



•ij, 



we obtain 



[d k±^d kj_g\A jg/pt{Xg,k±g)j^/f,[x^,kj 



dydM^d^q^ dydM^d^QT s J '^ '"^'"'^ '""^^"^ ''^'" '"''' ■-^^^7/ev-7, -±7J 

X 6\k^, + k^g - q^)a^,'^''%M') (5) 



and 



d^a^ d^a^ 

+ 



— [d k±^d k±g\jg/p[Xg, k±g)jy/f,[Xj, k_ 



dydM-^d^q^ dydM^d^qj. 

X 5\k^, + k^g - q^)al'^'\M') (6) 



where 



and the partonic cross section is given by 24 1 



1 2 47ra;asr/ 1 2\ i 1 + V^" 



V 



aV'^'^'iM') = i^el^^[{l + V - -v') In ^ _^^^—^ " (1 + v)VT^]. (8) 

Jiere, v = —— and M = s. 

Integrating eqs. dH) and (l6|) over M^, we obtain the difference and sum of - — - — and 

dyd^qrp 

- — - — for J /if) production. 
dyd^qrp 

Sivers asymmetry is defined as 

.sin{<>,y-0s) ^ / d(pq^[da^ - da^] sin(0g^ - (ps) , ■. 

^ j d(Pq^[da'^ + da^] ^' 

where da^ is differential cross section in qx or y variable and 0^^ and (f)s are the azimuthal 
angles of the J /if) and proton spin respectively. The weight factor in the numerator projects 
out the Sivers asymmetry. To evaluate asymmetry in y-distribution, we substitute 

da^ - da^= [ d(Pg^ f QTdqT [ "^ [dM^] [ [d^k^g]A'' fg/pr{xg,k^g) 



4mJ; 



X 



f-y/eiXy, Qt - f^±g) CToiM"^) sin(0,^ - (f)s) (10) 



and 



da^ + da^ = 2 d<i)q^ / gr dqr / \dM^] / [d'^k^g\fg/p{xg, k^g) 

J J Jiml J 

X f^/e{x^, Qt - k±g) cro{M^). (11) 

Thus at LO, the SSA depends on Weizsacker- Williams function, gluon distribution func- 
tion and gluon Sivers function. 

k I n dependence of the unpolarized pdf's is given by a simple factorized and Gaussian 
form [25]: 

fg/p{Xg, k^g) = h/pi^9)-7^^e-''^^/^'"-^\ (12) 

Similarly, transverse momentum dependent WW function is given by : 

f^/eix^. k^^) = f^/e{x.^^T^r^e-^'^'</^^'^'<\ (13) 



In Ref.[ll|, we have used the Sivers function given by|26|. 



where the gluon Sivers function, A'^ fg/p^-{xg) is defined as 

^^fg/pt{Xg) = 2Afg{Xg) fg/p{Xg). 

Here, J\fg{xg) is an x-dependent normahzation for gluon and 



COS( 0fc 



(14) 



(15) 



hik 



-Lo 



Ml 



(16) 



Ml is parameter obtained by fitting the recent experimental data corresponding to pion and 
kaon production at HERMES and COMPASS. 

The parameterizations for quark Sivers function Mu{x) and Md{x) have been fitted from 
SIDIS data and are given by [27|, 



Mfix) = NfX^'fn - xp ^^ ^\ . 



(17) 



Here af,bf, Nf for u and d quarks are free parameters obtained by fitting the data. However, 
ihere is no information available on J\fg (x) . In our analysis we have used two parametrizations 

27| 



(a) Afg{x) = iAfu{x)+Afdix))/2 

(b) Afg{x)=Afa{x). 



B. SCALE EVOLUTION OF TMDs 



In our previous analysis, we had chosen the scale of PDFs and Sivers function to be s. The 
TMD PDFs and Sivers function were obtained through DGLAP evolution by considering the 
evolution of factorized collinear part. The issue of scale dependence of TMDs has attracted 
a lot of attention in recent past and a new formalism for Q^ -evolution of TMDs has been 



developed 



|l6. 
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28| . Based on this TMD evolution, best fits of the SIDIS Sivers asymmetry 



have been performed and compared with earlier estimates extracted without using TMD 



evolution 



22|. In this work, we employ the strategy used in Ref. [22|] to take into account 



the TMD evolution of PDFs and the Sivers function and compare the resulting asymmetries 



in J /ip production with our earlier estimates which were obtained using DGLAP evolution. 
Since we are following the formalism of Ref . 22| , we will be brief in our description of the 



same. The details can be found in Ref. 



22|. 



In this formalism, Q^ evolution of /c^ -dependent distribution function is given by 22 1 



fqip{x, k±] Q) = fg/p{x, Qo) R{Q, Qo 



-kl/w^ 



TTW^ 



(18) 



where fq/p(x, Qo) is the usual integrated PDF evaluated at the initial scale Qo and w^ = 
w'^{Q,Qo) is the "evolving" Gaussian width, defined as: 

Q 



w\Q,Qo) = {ki) + 2g2\n 



Qo 



(19) 



Here, the evolution factor R{Q, Qo) is the limiting value of a function R(0, Qo, &t) that 
drives the Q^-evolution of TMDs in coordinate space and is given by 17|, l22| 



RiQ, Qo, br) ^ exp (in ^ /'' ^7i.(/i') + /'' %f ( fi, %) ] (20) 



where bx is the parton impact parameter and 



fJ-b 



b*{bj 



b*{b'2 



1 + b%lb 



2 Ih2 ■ 
max 



(21) 



with Ci = 2e-T^ and 7^ = -0.577 |29|. 

In the limit 67- — )■ 00 , R{Q, Qo, &t) -^ R{Q, Qo) and 6* -^ bmax- If and jk are anomalous 
dimensions which are given, at 0{as), by 



7f U^; —) = as[^^) — - - In — 
H^ vr \2 /i^ 



7/^(/i) = as(/i) 



2C. 



The TMD evolved Sivers function is given by 22 1 



TT 



,JV7 /^ L .r^\_ ^-^ ./^(^s) AN. 



^-kl/wl 



A^V,/pt(x,fc^;Q) = ^V2e^A^V,/pt(a;,Qo)i?(Q,Qo) ^^4 ' 



with 



where 



wl{Q,Qo) = {kl) + 2g2\n 



1 1 1 

+ 



Q_ 

Qo 



{kl) Ml {kl) ■ 



(22) 
(23) 

(24) 

(25) 
(26) 



{k\g) = 0.25 GeV'^, ^2 = 0.68, b^ax = 0.5 GeV-\ (27) 



Here, Mi is a best fit parameter |22j, Q^ = _§ and Qg ^ 1.0 GeV^ [29 1. 



C. ASYMMETRY IN J/i; PRODUCTION USING EVOLVED TMDs 

Using the TMD evolved PDF and Sivers function given in eq. (fTS!) and eq. i^^, and 
following the procedure in Ref.ful], we obtain the expression for numerator of asymmetry as 



dydPqrp dyd?qj. s 






((fc/))=exp[-#/(u.f + {fci,))l 



^'^<'^' '^"^ "'^"^H + ^YK)'^'" ^°^<^"' *"" <''^' ''" 



29| and expression for denominator as 



4mr 
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X i?(Q, Qo) ^^^2^(fc^ )] < (^ )• (29) 

We will use these expressions to estimate the asymmetry in the next section. It is known 
the Gaussian approximation is suitable only for low k± values as the Sivers function develops 



a tail at high fc^ values |l 81]. Therefore, while integrating over A;^, one needs to set the upper 
limit of integration below the value where the Sivers function starts deviating from the 
gaussian form. However, we are justified in integrating over the whole range of fc^ of gluon 
as the finite width of Gaussian dependence of WW function, in fact, provides additional 
damping, making the effective range of integration to be less than the value at which the 
above mentioned tail starts. We have, in fact, also performed the k^g integration in eq. 
f lTO|) and in eq. f lTTj) numerically over a finite range allowed by Gaussian approximation and 
verified that the result is the same as obtained by analytical integration. 

III. NUMERICAL ESTIMATES FOR THE ASYMMETRY IN J/V' PRODUCTION 
USING TMD EVOLVED SIVERS FUNCTION 

We will now estimate the magnitude of asymmetry using TMD evolved PDFs and Sivers 
function and compare the results with our previous estimates of asymmetry using DGLAP 

8 



evolution ll[|. We will also compare our results with asymmetry calculated using DGLAP 
evolution with another set of parameter extracted from DGLAP fits at Qq = 1.0 GeV. 
For the DGLAP evolution, we have estimated asymmetry using two different sets of 



parameters. The first set, w^ 
which we have used in Ref. 



lich we call DGLAP 1, consists of the values of best fit parameters 



ll|: 



Nu = 0.40, a„ = 0.35, 6„ = 2.6 , 

Nd = -0.97, aa = 0.44, ba = 0.90 , 

M^ = 0.19 GeV^. 



(30) 



These parameters are from new HERMES and COMPASS data 30|, l3l| fitted at Q^ 
2.4 Gel^^y _ 

The second set of parameters, which we call DGLAP2, has been extracted from DGLAP 
fits at Qo = 1.0 GeV 



22(: 



Nu = 0.45, an = 1.08, bu = 6.9 , 

Nd = -1.00, a^ = 1.7, bd = 6.9 , 

M^ = 0.19 Gel/l 



(31) 



In both cases, the value of {k'^g) is chosen to be the same as {k\) for quarks obtained 



in Ref. [32| by analysis of Cahn effect in unpolarized SIDIS from data collected in different 
energy and Q^ ranges assuming a constant Gaussian width. The value of {kj_^) has been 
chosen to be comparable to {k\). 

For TMD evolved Sivers function, we have used the parameter set fitted at Qo = 1 GeV 



given in Ref. 



21. 



a-u 



Nu = 0.75, Nd = -1.00, 6 = 4.0 ; 
0.82, ttd = 1.36, M^ = 0.34 GeV'^ 



(32) 



We have estimated the asymmetry with both kinds of parameterizations (labeled (a) 



and (b)). The estimates are obtained using GRV98LO for unpolarized 
functions |33| and Weizsaker- Williams function for photon distribution |34l ]. 
In figures 1-6 we have performed comparative study of Sivers asymmetry A 



uon distribution 



sin((/),j,-0s) 

N 



for 



DGLAP evolution and TMD evolution as a function of rapidity y and qt, the transverse 
momentum of J/^jj respectively for JLab (v^ = 4.7 GeV), HERMES (^i = 7.2 GeV), 



COMPASS (v^ = 17.33 GeV) and eRHIC (^i = 31.6 GeV and ^/^ = 158.1 GeV) energies. 
Before we start with a detailed discussion, we make one observation. In all the cases the 
asymmetries are smaller than our earlier estimates. However, they still remain sizable, 
except for the lowest energy at JLab, where any way it is a difficult measurement due to the 
low event rate. 

In figure 1-5, SSA with DGLAP evolution with parameter sets DGLAPl and DGLAP2 
and TMD evolution is compared using parametrisation (b). To get the qx distribution of 
asymmetry we have integrated over all possible values of y and to get the y distribution of 
asymmetry we have integrated over qt from to 1.0 GeV. The masses of the c quark and 
the D meson are taken to be 1.275 and 1.864 GeV, respectively [_35|. In TMD evolved PDFs 
and Sivers function, the initial scale Q?. has been chosen to be 1.0 GeV as we are using Sivers 
funct.on pa.a,„ete. of Ref. ^ fl«ed at th. scale. Since we a.e u.„g the cole, evaporation 
model, evolved TMDs are evaluated at Q^ = s which varies between 4ml ^^"^ Amj^ and is 
the only relevant scale for J/ip production in the GEM. 

We note that all the asymmetries are suppressed with respect to our earlier predictions 
using the older model and parameter set DGLAPl ll|]. To be specific, we see that for pa- 
rameterization (b) of gluon Sivers function, with TMD evolution, the maximum asymmetry 
in y distribution is reduced from approximately 15% to 5% at different rapidity values, for 
energies of JLab, HERMES, GOMPASS, eRHIG-1 and eRHIG-2 experiments. For qx dis- 
tribution of asymmetry with parameterisation (b), maximum asymmetry is reduced from 
approximately 15% to 5% at gr=0.75 GeV for energies of JLab, HERMES, GOMPASS and 
eRHIG-1 experiments. For higher eRHIG energy, asymmetry in g^ distribution is reduced 
from approximately 5% to 1%. 

We have also compared our estimates of asymmetry using TMD evolved PDF's with 
asymmetry calculated using DGLAP evolved PDF's with parameter set DGLAP2. The 
parameters in these two cases have been extracted from TMD fits at Qq = 1.0 GeV and 
DGLAP fits at Qq = 1.0 GeV respectively assuming a Gaussian width of 0.5 GeV at this 
initial scale. In this case also, we find appreciable suppression but the difference between 
DGLAP and TMD estimates is smaller now. In particular, the asymmetry is not changing 
much for the TMD evolution and DGLAP evolution in case of y distribution of JLab exper- 
iments. For qt distribution of asymmetry with parametrisation (b), maximum asymmetry 
is reduced from approximately 8% to 5% for eRHIG-1 experiment. For qx distribution of 

10 



asymmetry in eRHIC-2 experiment, asymmetry is reduced from approximately 1.5% to 1%. 
In general, we find that the estimates using this parameter set are closer to the estimates 
using TMD evolved PDF's as compared to our earlier estimates using parameter set fitted at 
Q^ = 2.4 GeV^ . This seems to suggest the importance of Q^ dependence of Gaussian width 
as all the earlier DGLAP fits of Sivers function have been obtained assuming a constant, 
scale-independent Gaussian width. 

In figure 6, we have compared the SSA for parameterization (a) and (b) using TMD evolu- 
tion. As in case of DGLAP evolution, we find that the asymmetry is higher for parametriza- 
tion (b) than that for parametrization (a) in case of TMD evolution as well. 

IV. DISCUSSION AND SUMMARY 

In this work, we have given numerical estimates for J/ip production in low virtuality 
electro-production at JLAB, COMPASS, HERMES and eRHIC energies, using the TMD 
evolved parton distribution functions and Sivers function. At leading order, this asymmetry 
is a clean probe of the gluon Sivers function. We use here the analytic formulation for 
incorporating the scale evolution of the TMDs and calculate the asymmetries in the color 
evaporation model for J/ip production. We compared these estimates with those obtained 
by us previously using a different parameterisation wherein the Q^ dependence came from 
DGLAP evolution of the unpolarised gluon densities. The present estimates, while reduced 
substantially, when compared to our earlier ones, are still sizable for most experiments. 

Noting that the average Q^ for our process is ~ 10 GeV^, the suppression we find is seen 
consistent with the suppression in going from HERMES to COMPASS energies, obtained 
in Ref. 19| in the context of SIDIS. It should be also remembered that in case of J/^ 
production, the transverse momentum dependence of the WW function, can also have non 
trivial effects. 

In the model under consideration, the amount of asymmetry is affected by the evolution 
factor R{Q, Qq) as well as by the Q^-dependence of Gaussian width. In color evaporation 
model, Q^ varies between 4?n^ and Am'jy and R{Q, Qo) does not vary much over this range. 
Hence, this factor approximately gets canceled between the numerator and denominator 
in the evaluation of asymmetry. Thus, the noticeable suppression comes from logarithmic 
dependence on Q^ of the Gaussian width. This is consistent with the observation in 22[ that 

11 



the substantial decrease in asymmetry is mainly due to Q^ dependence of the Gaussian width 
of TMD PDFs. Hence it is important to understand this Q^ - dependence of Gaussian width 
in order to be able to make predictions which can be compared with the future experiments. 
In case of J /ip production, there is yet another convolution with transverse momentum 
dependent WW function. In principle, one needs to model its scale dependence as well. In 
the end it should be noted that, since the hard scale is the same for all beam energies, the 
Sivers function is probed in the same Q^ range. This explains similar amount of reduction 
of asymmetries seen at all energies. 

Our estimates in the present work are based on the approximate form of Sivers function in 
Eq. (!24l) which was obtained in Ref. [2^ using TMD evolution equations of the Refs. [l6|-[l8|. 



This form, called the analytic form in Ref. 2^, is obtained by making an approximation 



6t — !■ oo in Eq. (120!) . An exact form can be obtained from Eq. (!20l) by first evaluating the 
TMD evolved PDFs in coordinate space -F(x, hx'-, Q) and then taking their Fourier transform. 
However, the approximate form we have used is valid for 67- > l.QGeV~^ if the typical k± 
involved is of the order of 1 GeV. Using the exact form may lead to slight differences from 
the present estimates and we plan to address this issue in future work. 
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JLab, VF=4.7 GeV 
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FIG. 1. The Sivers asymmetry ^'^^^''t ^s) ^^^ ^^^t ^ e+J/il)+X at JLab energy (^i = 4.7 GeV) 
as a function of y (left panel) and qt (right panel) for parametrization (b). The solid (blue) line 
corresponds to results obtained using TMD evolution. The dashed (red) and dotted (black) line 
corresponds to DGLAP evolution with DGLAP fit parameters at Qq = \J'1A GeV and Qq = 1 GeV 
respectively . The integration ranges are (0 < gy < 1) GeV and (—0.25 <y< 0.25). 
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FIG. 2. The Sivers asymmetry A^'^^^"'^ ^^^ iov e + p'^ ^ e + J/xp + X at HEMRES energy 
(-y/s = 7.2 GeV) as a function of y (left panel) and qt (right panel) for parametrization (b). The 
solid (blue) line corresponds to results obtained using TMD evolution. The dashed (red) and dotted 
(black) line corresponds to DGLAP evolution with DGLAP fit parameters at Qq = \pTA GeV and 
Qo = 1 GeV respectively . The integration ranges are (0 < qr < 1) GeV and (—0.6 < y < 0.6). 
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COMPASS, Vs"=17.33 GeV 
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FIG. 3. The Sivers asymmetry ^^"^^"t -^^v f^j. g _^ ^t _^ g + J/ip + X at COMPASS energy 
{^/s = 17.33 GeV) as a function of y (left panel) and qt (right panel) for parametrization (b). The 
solid (blue) line corresponds to results obtained using TMD evolution. The dashed (red) and dotted 
(black) line corresponds to DGLAP evolution with DGLAP fit parameters at Qq = \pTA GeV and 
Qo = 1 GeV respectively . The integration ranges are (0 < g^^ < 1) GeV and (—1.5 <y< 1.5). 
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FIG. 4. The Sivers asymmetry aJ"^*^^ -^s) ^^ g ^ ^t ^ g + J/^ + X at eRHIC-1 energy (^ = 
31.6 GeV) as a function of y (left panel) and qr (right panel) for parametrization (b). The solid 
(blue) line corresponds to results obtained using TMD evolution. The dashed (red) and dotted 
(black) line corresponds to DGLAP evolution with DGLAP fit parameters at Qq = V2.4 GeV and 
Qo = 1 GeV respectively . The integration ranges are (0 < qr < 1) GeV and (—2.1 < y < 2.1). 
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eRHIC-2, Vr=158.1GeV 



eRHIC-2, Vs"=158.1 GeV 
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FIG. 5. The Sivers asymmetry ^^'^('^^t-'^s) ^^ g ^^t ^ g + J/?/' + X at eRHIC-2 energy (^ = 
158.1 GeV) as a function of y (left panel) and qt (right panel) for parametrization (b). The solid 
(blue) line corresponds to results obtained using TMD evolution. The dashed (red) and dotted 
(black) line corresponds to DGLAP evolution with DGLAP fit parameters at Qq = \pTA GeV and 
Qo = 1 GeV respectively . The integration ranges are (0 < g^ < 1) GeV and (—3.7 < y < 3.7). 
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FIG. 6. The Sivers asymmetry ^^"^"^"t v^^v f^j. g _^ ^t _^ g + J/ip + X at COMPASS energy 
(-y/s = 17.33 GeV) as a function of y (left panel) and qt (right panel) for TMD evolution. The 
solid (red) line corresponds to results obtained using parametrization (a) and the dashed (blue) 
line correspond to parametrization (b). The integration ranges are (0 < gr < 1) GeV and (—1.5 < 
y<1.5). 
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